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Abstract 

The treatment of degenerate states within Kohn-Sham density functional theory (KS-DFT) is 
a problem of longstanding and current interest. We propose a solution to this mapping from the 
interacting degenerate system to that of the noninteracting fermion model whereby the equivalent 
density and energy are obtained via the unifying physical framework of quantal density functional 
theory (Q-DFT). We describe the Q-DFT of both ground and excited degenerate states, and for the 
cases of both pure state and ensemble v-representable densities. The Q-DFT description further 
provides a rigorous physical interpretation of the corresponding KS-DFT energy functionals of the 
density, ensemble density, bidensity and ensemble bidensity , and of their respective functional 
derivatives. We conclude with examples of the mappings within Q-DFT. 
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The treatment of degenerate states within the context of Kohn-Sham density functional 
theory KS-DFT [l] is a problem of longstanding [2] and continued recent (3, |^| interest. 
The basic idea underlying the KS-DFT methodology is the mapping from the interacting 
electronic system as described by Schrodinger theory to that of a system of noninteracting 
fermions such that the equivalent density and total energy are obtained. The existence of 
the model system is an assumption. Further, there is the distinction between the S system 
of noninteracting fermions whereby the equivalent density is obtained from a single Slater 
determinant, and the noninteracting system whose orbitals could be degenerate so that 
the density is obtained from a weighted sum of the Slater determinants constructed from 
the orbitals. The following cases have been considered. The mapping from: (a) a pure 
degenerate ground state |5j; (b) a pure degenerate excited state |2J. In addition, maps to 
obtain the density and energy constructed from (c) an ensemble of pure degenerate ground 

n 

statesfj], and (d) an ensemble of pure degenerate excited states[3|, have been developed. 
The interest in the ensemble cases stems from a ground state theorem due to Levy [|| 
and Lieb jl). According to the theorem, most ensemble densities constructed from pure 
degenerate ground states are not interacting v-representable. In other words, no single 
ground state wave function of the Schrodinger Hamiltonian will yield this ensemble density. 
Such ensemble densities are said to be ensemble v-representable. The translation of the 
theorem to the S system [f| means that there is no single Slater determinant that leads 
to this ensemble density. At this time, the question of interacting v-representability of 
the ensemble density of degenerate excited states is still unanswered jsj]. In this paper we 
propose a solution to the problem of mapping from the interacting degenerate system to 
that of the equivalent noninteracting fermion model within the unifying physical framework 
of quantal density functional theory (Q-DFT). 



In contrast to the present work, the KS-DFT description of the noninteracting system 
is in terms of an energy functional of the density, and its functional derivative. For case 
(a), the energy is a functional of the degenerate ground state density; for (b), the energy 
is a bidensity functional of the ground and excited state densities, with the functional 
derivative taken at the excited state density; for (c), the energy, which is a functional of the 
ground state ensemble density, is constructed by the ensemble generalization of the coupling 
constant scheme; and for (d) the energy is also a bidensity functional, in this instance of 
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the ground state and excited state ensemble densities, with the functional derivative taken 
at the ensemble density. 

We describe here the Q-DFT of both ground and excited degenerate states, and for the 
cases of both t he p ure state and ensemble v-representable densities. Q-DFT is an alternate 
description 0, Q] of the mapping from the interacting to the noninteracting system. A 
key point to note is that within Q-DFT the densities concerned are always obtained from 
solution of the Schrodinger equation. Thus, the pure state density, and each component 
of the ensemble density, are interacting v-representable. The assumption of existence of 
an S system in Q-DFT therefore means that the pure state density and each component 
of the ensemble density are also noninteracting v-representable. (By an S system we mean 
a noninteracting fermion system whose wavefunction is a single Slater determinant, and 
which maybe in a ground or excited state). In Q-DFT, the local electron interaction 
potential energy v ee (r) of the S system, and the expression for the total energy are in 
terms of 'classical' fields whose quantal sources are quantum-mechanical expectations of 
Hermitian operators. We begin by (i) describing the Q-DFT of the individual degenerate 
pure state. For the mapping from a degenerate ground state of the interacting system, the 
corresponding S system is in its ground state. For the mapping from a degenerate excited 
state, the state of the S system is arbitrary in that it may be in a ground or excited state. 
In either case, the highest occupied eigenvalue is the negative of the ionization potential. 

n 

(The mapping from a nondegenerate excited state [10] is similar.) For the ground and 
excited state ensemble cases, we describe two different schemes within Q-DFT. Thus, 
(ii) in the first, the corresponding noninteracting system ensemble density is obtained by 
constructing g S systems, where g is the degeneracy of the state. Once again for excited 
states, the g S systems may either be in a ground or excited state or a combination of the 
two. Next, (Hi) we describe the Q-DFT whereby the ensemble density is obtained from a 
single noninteracting fermion system whose orbitals could be degenerate. The construction 
of this model system is a consequence of the linearity of the differential virial theorem. 
Here the highest occupied eigenvalue is degenerate, and the ensemble density is obtained 
from the resulting Slater determinants as described by Ullrich and Kohn [4] whose work in 
turn is based on that of Chayes et al |rjj. Again for the mapping from an excited state, 
the noninteracting system may be in a ground or excited state. The Q-DFT description 
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then (iv) provides the physics underlying all the various KS-DFT degenerate state energy 
density and bidensity functionals, and of their functional derivatives. Finally, (v) we 
present examples that demonstrate the above mappings within Q-DFT. 

A. The Q-DFT of the bound individual degenerate pure state is as follows. The 
Schrodinger equation for a degenerate state whether ground or excited is 

iftt n ,„(X) = [f + V + E/]tt n ,„(X) = £ n tt n ,„(X), (1) 

f = -|EiVj 2 , V = T,iv(ri), U = \ ]^~\, where *„ j7 ,(X) and E n are a bound 
degenerate state wave function and energy, n corresponds to the state, and r\ = 1, ...g n the 
degeneracy, X = x 1; ...x N ,x = rer, with a the spin coordinate. As the equations to follow 
are valid for arbitrary states, we drop the subscript n. The degenerate pure state density 
p v (r) = | p | *,,), where p = £i 8(r - r;), and the energy E v = \ H \ 

The corresponding differential equation for the S system of noninteracting fermions with 
the same density is 

[__V 2 + ^(r)+^(r)]^(x) = £i&(x); i = 1, -N, (2) 

with 

Pr ,(r) = ($,{0,} | p | $,{0,}) = ]T | 0,(x) | 2 , (3) 

i,cr 

and $^{0*} is the smg/e Slater determinant of the orbitals 0«(x). This is the S system wave 
function. The electron-interaction potential energy v ee ^{v) is representative of electron cor- 
relations due to the Pauli exclusion principle, Coulomb repulsion, and Correlation-Kinetic 
effects. Correlation-Kinetic contributions to the potential energy are a consequence of the 
difference in kinetic energy between the interacting and noninteracting systems. The poten- 
tial energy v ee>ri (r) is the work done to move a model Fermion in the force of a conservative 
field ^j(r): 

VeeM = - [ r ^(r')-dY, (4) 

J oo 

where ^(r) = S ee ^(r) + Z tctV (r). The fields £ eetV (r) and Z tc ^{r) are not necessarily conser- 
vative. Their sum always is. The electron-interaction field £ eejV (r) is representative of Pauli 
and Coulomb correlation: £ ee , ri (r) = e ee , ri (r)/ p ri (r), where the electron-interaction 'force' 
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e ee)r? (r) is obtained via Coulomb's law as e e6)r? (r) = / dr'P^(rr')(r — r')/ | r — r' | 3 , 
where P„(rr') = (^|P(rr')|^), with P(rr') = £y<J(r - r, ; )<5(r' - r^). Equiv- 
alently, the field £ ee)r? (r) may be thought of as being due to its quantal source, 
the pair-correlation density g v (rr') = P v (rr') / p n (r) . The Correlation-Kinetic field 
z t c , v (r) = Z SyV {r) - Z v (r), Z v (r) = z v (r; [j v ])/Pr,(r), Z s>v (r) = z SyV (r; ^J/p^r), 
and where Z v (r) and Z s ^(r) are the interacting and S system kinetic fields, respec- 
tively. The kinetic 'force' z r7 (r) is defined by its component z V}0 , = ZJ^pdta^/drp, 
with t a)1 g(r, [7^]) = f [0 2 /^ r a^ r |3 + d 2 / dr'pdr'^^r' , r")| r / =r // =r the kinetic energy density 
tensor. The source of the kinetic field Z v (r) is the spinless single particle density matrix 
7 ,(r,r') = (V V \X\%), X = A + iB, A= ± E#(r, - r)T,(a) + tffa - r')I}(-a)], B = 
— 2 Sj[^( r j ~~ r )^j( a ) ~~ ^( r j ~~ r ')^j(~ a )]; ^?'( a ) i s a translation operator, and a = r' — r. 
The field Z s (r) is defined in a similar manner in terms of the S system Dirac density matrix 
7 ^(r,r') = ($„{&}|A-|$„{&}) = ^M^rV). 

The proof of Eq.(4) follows by equating the differential virial theorems Q| for the 
interacting and S systems which are, respectively 

Vu(r) = -F„(r) and Vv(r) = -F,,„(r), (5) 

where F J? (r) = -5 ee ^(r)+P^(r)+2^(r), F Sjr) (r) = V^W+D^rj+^^r), the differential 
density field T> n (r) = d^(r)/p^(r), d r? (r) = — iVV 2 /?r/( r )- Thus, one obtains 

Vv ee>v (r) = -Fti(t), (6) 

from which the interpretation of Eq. (4) follows. 
The total energy of the degenerate state 77 is then 

E v = T S)V + J p v (r)v(r)dr + E ee>rj + T CtV , (7) 

where T s>v = {^> v {(p i }\T\^ rj {(j) i }) is the S system kinetic energy, and the electron-interaction 
E ee>r] and Correlation-Kinetic T c rj energies in terms of the respective fields are 

T c , v = ^Jdrp v (r)r-Z tc>v (r). (9) 
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These expressions are independent of whether the fields £ eeir? (r) and Z tc;V (r) are conservative 
or not. 

The S system whereby the density and total energy equivalent to that of the interacting 
system degenerate state 77 is defined by Eqs.(2)-(4) and (7). If the degenerate state is 
excited , the S system may be constructed to be either in a ground or excited state. Since 
the electron-interaction field remains unchanged, the difference between the corresponding 
potential energies is independent of the Pauli principle and Coulomb repulsion and due 
entirely to the corresponding Correlation-Kinetic fields Z tc , n (r). Hence, the potential energy 
v eetV (r) is different depending on whether the S system is in a ground or excited state. In 
either case, the highest occupied eigenvalue of the S system differential equation is the 
negative of the ionization potential. This follows by equating the asymptotic structure of 
the density for the interacting and S systems. 

In the transformation from an excited pure degenerate state to an S system in its ground 
state, the fact that the interacting system wave function has nodes is of no relevance. By con- 
struction, the S and interacting system density p(r) are equivalent, and the density p(r) > 0. 
Such a mapping for an excited pure nondegenerate state has been demonstrated in Ref. [10]. 

B. We next describe the first of two ways of obtaining the ensemble density and energy of 
the degenerate states via Q-DFT. The interacting system ensemble density matrix operator 
D(XX') is defined as 

J D(XX') = E^;( X )^(X / ); X>„ = 1; 0<^<1, (10) 

rj=l r/=l 

so that the ensemble density p e nsij) and energy E ens are respectively 

9 

Pensij) = tr(Dp) = J2 VnPv( T )> ( n ) 
v =i 

and 

9 

E ens = tr(DH) = J2 "v E vi (12) 
n=i 

with p v (r) and E n as defined previously. (There are ensemble densities that cannot be 
represented by a single Slater determinant. However, its pure state component density can 
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always be reproduced by an S system). 



For each degenerate state r], the density p r? (r) and energy E v can be constructed from 
an S system as described in part A. Thus, the ensemble density and energy of Eqs. (11) 
and (12) may be obtained from g S systems. Each S system contributing to the ensemble 
density may be in a ground or excited state. Note that the electron-interaction potential 
energy v ee ^(r) for each of the g S systems will be different. Further, v ee ^ v (r) will be different 
depending on whether the particular S system is in a ground or excited state as explained 
previously. Thus, the ensemble density and energy within Q-DFT are obtained by replacing 
the p v (r) and E v on the right hand sides of Eqs. (11) and (12) by the corresponding S 
system equivalents of Eqs. (3) and (7), respectively. 



C. The ensemble density and energy may also be determined from a noninteracting 
fermion system whose orbitals could be degenerate as constructed within Q-DFT. According 
to Chayes et al [ill, the ground state ensemble density may be determined as a unique 
weighted sum of squares of a finite number g of degenerate wave functions of this system. 
The potential energy v ee (r) of these noninteracting fermions is then determined via Q-DFT 
as follows. Rewrite the interacting and noninteracting system differential virial theorems of 
Eq. (5) as 

g 

p ens {r)Vv (r) = 5^w^(r), (13) 

77=1 

where f^(r) = e eefl (r) - d^(r) - z v (r), and 

9 

p ens (r)Vv(r) = ^2uJ v f StV (r), (14) 

,7=1 

where f Sl77 (r) = p n {r)Vv eetri {Y) — d^(r) — z S)7? (r). Equating Eqs. (13) and (14) leads to 

9 9 

Ur,Pr,(r)Vv ee>v (r) = w T?qjj(r), (15) 

r]=l r]=l 

where q^(r) = e eefl (r) + z tc;V (r), z tej7? (r) = z.J r) - z v (r). Eq. (15) is a consequence of 
the linearity of the differential virial theorem [12] . As we require a single effective potential 
energy v s (r) = v(r) + f ee (r), we replace v eeiV (r), in Eq. (15) by v ee (r) to obtain 

Vt; ee (r) = -Q(r), (16) 
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where Q(r) = —(J2f ] =i UJ rj ( irj(' r ))/Pens( r ), Thus, the electron-interaction potential energy 
v ee (r) is the work done in the conservative field Q(r): 

v ee (v) = - [ r Q(r')-dY. (17) 

Note that the components q ?? (r) are conservative so that Q(r) is conservative, and hence 
v ee (r) is path independent. 

For the occupation of orbitals we follow Ullrich-Kohn [4J. Accordingly, all levels are 
occupied except the highest (h) which is q- fold degenerate and partially occupied. The 
number of the model fermions in these levels are N h < 2q. The ensemble density which is 
a weighted sum of the degenerate Slater determinants is 

N-N h q 
Pens(r)= £ I^(X)| 2 + E/^"(^^)| 2 , (18) 

with < fi < 1, and fi = X^=i where 9i tV = 1 if the orbital 0^(x;i?) occurs in 

the determinant Q v {<j)i}, and otherwise. Here the <^(x; R) are appropriately rotated 
(R) orbitals determined self-consistently together with the fi and the lower lying orbitals 
leading to the ensemble density. The ensemble energy is obtained from the g Slater deter- 
minants as in part B. Once again for an excited state ensemble density, the corresponding 
noninteracting system may be in a ground or excited state. 

Note that the methodology of construction of the g S systems of part B also follows 
from Eq. (15). We believe that it is easier to construct the g S systems of part B than it 
is to construct the single noninteracting system of part C. This is because each of the g S 
systems may be constructed independently. 

From the above degenerate state Q-DFT description it is then possible to provide 
a rigorous physical interpretation for each energy functional and functional derivative 
of the corresponding KS-DFT. In each case, the local potential energy of the model 
fermions is the work done in a conservative field. The energy in turn may be expressed 
in terms of the components of this field. Thus, for example, the KS-DFT degenerate 
ground state electron-interaction energy function, E M of t he — de n - 
sity is the ensemble sum of the electron-interaction E ee>ri and Correlation-Kinetic T c?? 
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energies. The functional derivative is the work done to move the model fermion in the 
conservative field Q(r). The same interpretation applies to the bidensity energy func- 
tional E^ e s [p gr , p ens ] |3| of degenerate excited state KS-DFT, and of its functional derivative. 



ae demon- 
This 
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The Q-DFT mapping from a pure degenerate excited state to an S system can 
strated via the first excited triplet state of the exactly solvable Hooke's atom 
atom is comprised of two electrons with a harmonic external potential energy v(r) = \ur 2 . 
The triplet state wave functions are of the form 



*( ri r 2 ) = a o e-^%-- 2 /1l + C i y|r + C 2 (|)r 2 + C'3(|) 3 / 2 r 3 ]nn l ^0), (19) 

where Y^ m (8,(f)) is the spherical harmonic with I = l,m = — 1,0, l,r = r 2 — r 1; R = 
(r 2 + i"i)/2, and Co, Ci, C 2 , C3 are constants. For each degenerate wave function, the 



corresponding S system, in either a ground or excited state, can 



2e_ determined in a manner 
These results will be 



10]. 



similar to the transformation of the first excited singlet state 
presented elsewhere. The ensemble density of these degenerate states then follows from the 
3 S systems. (Note that the ensemble density for this two electron model is v-representable. 
However, the methodology for constructing the g S systems is the same whether or not 
the density is v-representable.) Another example is that of the noninteracting Be atom j^J. 
Here the ensemble density, which is not v-representable, is the weighted sum of the density 
of 4 S systems in the states ls 2 2s 2 , ls 2 2pj (i = x,y,x). This latter model atom is also an 
example 4] of the noninteracting fermion system that leads to the ensemble density with 
appropriately rotated highest occupied orbitals. 



In conclusion, we have described via Q-DFT the physics of mapping from a degenerate 
state of Schrodinger theory to that of a model system of noninteracting fermions such that 
the equivalent density and energy are determined. The cases of both pure and ensemble 
v-representable densities are explained. The framework is general and formally the same 
for both degenerate ground and excited states. 
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